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Nucleation of colloids and macromolecules: does the nucleation pathway 
matter? 
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A recent description of difFusion-limited nucleation based on fluctuating hydrodynamics that extends classical 
nucleation theory predicts a very non-classical two-step scenario whereby nucleation is most likely to occur 
in spatially-extended, low-amplitude density fluctuations. In this paper, it is shown how the formalism can 
be used to determine the maximum probability of observing any proposed nucleation pathway, thus allowing 
one to address the question as to their relative likelihood, including of the newly proposed pathway compared 
to classical scenarios. Calculations are presented for the nucleation of high-concentration bubbles in a low- 
concentration solution of globular proteins and it is found that the relative probabilities (new theory compared 
to classical result) for reaching a critical nucleus containing Nc molecules scales as e"^"/'^ thus indicating 
that for all but the smallest nuclei, the classical scenario is extremely unlikely. 



I. INTRODUCTION 

Nucleation - whether homogeneous or heterogeneous 
- is a paradigmatic example of self assembly. It occurs 
when a physical system can be in two or more states 
that are separated by free-energy barriers. If the system 
is initially in one of the metastable states, thermal fluc- 
tuations can drive it over the free energy barrier and into 
a more stable state. Because the energy needed to over- 
come the barrier scales with the spatial size of the system, 
the process occurs locally via the formation of a finite 
sized cluster or nucleus. In many processes of interest 
such as the crystallization of proteins from solutionii^, 
the formation of snowflakes^ and the crystallization of 
polymorphic solids^ the system may pass through one or 
more intermediate metastable states before arriving at 
the final, stable state. In fact, the heuristic known as 
Ostwald's rule of stages specifically states that a system 
will pass in turn from one state to another having the 
next lowest free energy until it reaches the minimal en- 
ergy state. Various arguments can be given in support 
of this rule^ and they can be grouped into two classes: 
either it arises because the barriers separated "similar" 
states are smaller than those separating disparate states 
or it arises due to the kinetics of the transition. The for- 
mer reason is an application of another heuristic known 
as the Stranski-Totomanow conjecture which states that 
the observed transition will be the one corresponding to 
the minimal energy barrier—. The latter is more difficult 
to characterize as it may depend on microscopic details 
of the various states and the dynamics of the system. 
Recent work on model systems indicates that, depending 
on the free energy landscape and the dynamics, either 
Ostwald or Stranski-Totomanow or both may be correct 
or noti. 

The common issue at question is the description of the 
nucleation pathway. Clearly, a theoretical description for 
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the nucleation pathway that goes beyond the empirical 
heuristics mentioned above must be based on a dynam- 
ical description of the nucleation including the role of 
thermal fiuctuations, mass and energy transport and the 
structure of the various phases. A framework for such 
a dynamical formulation of nucleation has recently been 
described^ii^. For large clusters, it was shown to repro- 
duce classical nucleation theory (CNT) in the weak-noise 
limit. For a diffusion-limited dynamics - appropriate for 
the description of colloids and macromolecules in solution 
- it was shown that the relative probability of different nu- 
cleation pathways could be easily calculated and that the 
most likely path (MLP) could be determined by steepest- 
descent on the free energy surface. This framework can 
be contrasted with numerous proposals for determining 
the nucleation pathway based solely on the properties 
of the free energy surface with no dynamical input, see 
e.g. [lol - [l6l . As discussed below, because of the heuristic 
nature of the latter, evaluating their relative merits has 
proven difficult. 

The purpose of this paper is to illustrate the compari- 
son of different candidate nucleation pathways by means 
of their relative probability. The particular example in- 
vestigated is the formation of high concentration droplets 
in a low-concentration solution of globular proteins. This 
process is analogous to the vapor-liquid transition in sim- 
ple fluids and, while being of intrinsic interest as part of 
the process of crystallization in globular proteins, has 
the practical advantage of allowing for a relatively sim- 
ple theoretical description as described in Ref. d and d 
and below. For this problem, the recently developed dy- 
namical theory of nucleation predicts a very different nu- 
cleation pathway than do the older DFT-based theories. 
The new prediction is that the process of droplet nucle- 
ation involves two steps: first, a long- wavelength density 
fiuctuation forms and then a nucleation event takes place 
within this fiuctuation. The older, classical view is that 
nucleation begins with a spatially localized cluster that 
grows monotonically. From a point of view focused on 
the comparison of the theories, it is therefore of interest 
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to compare the relative likelihood of these two pathways 
so as to determine whether the difference between them 
is qualitative or quantitative. More broadly, in the con- 
text of the general theory of nucleation this allows us to 
address the question of how important it is to choose the 
"right" path when there are multiple possible paths that 
pass through the (unique) critical cluster. Note that a 
feature of the heuristic rules is that they are phrased en- 
tirely in terms of the free energy maxima (the barriers) 
and minima (the states) and, hence, the provide no guid- 
ance on this question. In this sense, they are in accord 
with the common intuition that all that matters in nu- 
cleation are the free-energy extrema which, in addition 
to their use in choosing pathways, are also the only rele- 
vant quantities entering the Classical Nucleation Theory 
(CNT) for nucleation ratesii. Here, it will be shown 
that even paths beginning and ending on the same states 
and passing over the same barriers can have wildly dif- 
ferent probabilities of occurrence. The paths compared 
will be the most likely path as determined from the dy- 
namical theory and a path determined from one of the 
non-dynamical, free-energy methods which, by means of 
direct comparison to simulationi^, is known to give a very 
accurate description of the free energy barrier for nucle- 
ation in proteins. Since this method has also recently 
been applied to the study of nucleation of wetting films 
on curved substratesi^ and of ordered phases of block 
copolymersi^, where multiple candidate pathways were 
found, it is of particular interest to use as a test-case. 
This paper therefore serves two puposes: (i) to introduce 
a method of comparing the likelihood of candidate nucle- 
ation pathways, however they are arrived at; and (ii)to 
use this method to determine whether the proposed non- 
classical pathway is significantly different (more proba- 
ble) than a "classical" alternative. 

In Section II, the elements of the theoretical descrip- 
tion are reviewed. Section III describes the application 
to the nucleation of protein-rich droplets in solution, in- 
cluding a detailed comparison of the different candidate 
pathways and a computation of their relative probabili- 
ties. The paper ends with a brief discussion of our con- 
clusions. 



son, the concentration of the large particles is equivalent 
to their density and the two terms will be used inter- 
changeable in the following. Further details of the micro- 
scopic model can be found in Ref. [^. The fundamental 
quantity with which the phase transition is character- 
ized is then the local concentration (or number density) , 
/9(r). The density /concentration is commonly assumed 
to be spherically symmetric and this assumption will be 
used throughout the present development. The initial, 
metastable, system is characterized by a uniform density 
p (r) = Pi where pi is a minimum of the bulk free energy. 
The new phase is also characterized by a uniform den- 
sity, p{r) = pf, where p/ is the global minimum of the 
bulk free energy. We assume throughout the existence of 
a Helmholtz free energy functional, F [p], so that, e.g., in 
the grand canonical ensemble the appropriate free energy 
is ft[p] ~ F [p] — /iiV where = J p{r) dr is the total 
number of particles. For a uniform density the free energy 
becomes an ordinary function, F[p{r) = Pi] = F{pi), 
and the conditions for the uniform phases to be minima 
are the usual relations yF' (pi) = fJ, = yF' {pf). 

Nucleation proceeds by the formation of a cluster con- 
sisting of the new phase which grows until it consumes 
the entire system. Despite the fact that the new phase is 
energetically favored, small clusters are unstable due to 
the dominance of surface tension effects. When a clus- 
ter is sufficiently large, the lowering of the cluster energy 
due to increasing the size of the bulk region inside the 
cluster outweighs the cost of increasing the surface area 
and growth is favored. These regimes are separated by a 
saddle point in the free energy called the critical cluster 
which necessarily satisfies the relation 
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There are generally two approaches to representing the 
density function. One is to simply discretize space by 
introducing a lattice of points = iA so that one works 
with a series of values pi = p{iA). An alternative is to 
use parametrized functional form such as a hyperbolic 
tangent, 



II. THEORY 

The present theoretical development is concerns a col- 
lection of particles - molecules or colloidal particles - that 
interact with one another via a prescribed pair potential 
and which are also subject to random, Brownian forces. 
This is a simple model for large particles in a bath of 
small particles wherein the effect of the bath or solution 
is incorporated via the effective interaction between the 
large particles and random (Brownian) forces acting on 
the large particles. As such, the bath is not explicitly 
represented except through the amplitude of the random 
noise which in turn determines the (low-concentration) 
diffusion constant for the large particles. For this rca- 
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or an exponential form 
p{r) = 
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where, in both cases, there are four parameters: a radius, 
R, a width, w, the interior density, po that characterizes 
the central density for large {R » w) clusters and the 
density far from the cluster, p^a- It is easy to show that 
in the thermodynamic limit the latter must be a min- 
imum of the free energy so that we will normally have 
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that it is equal to the initial density = pi . The advan- 
tage of the first method, discretization, is that it is clear 
that one can approach the continuum limit by decreasing 
the lattice spacing whereas the advantage of the second 
method, parametrization, is that one can hope to get 
good results with relatively few parameters. In fact, the 
two methods can be viewed as two different approaches 
to parametrization and other possibilities - such as rep- 
resentation in terms of Fourier components or projection 
onto some other set of basis functions - have the same 
characteristic. We can therefore without loss of general- 
ity assume that the density field is represented by a col- 
lection of parameters denoted generically as Xt so that 
p{r) = f (r';x), for some function /. The change in the 
density profile as a function of time therefore becomes a 
change in the parameters so that more generally we have 
p (r; t) ^ f (r; x (t)). Thus, specification of the evolution 
of the parameters, x(t), corresponds to the prescription 
of a path in density space. The same notion holds when 
the path is parametrized by some other quantity rather 
than time - e.g. the equivalent of a reaction coordinate. 

To give a concrete illustration, Ghosh and Ghosh use 
the exponential profile with the radius acting as the reac- 
tion coordinate. They determine the other parameters by 
minimizing the free energy while holding the radius con- 
stant thus parametrize the path by the radius. However, 
this is not a unique prescription since one could equally 
well parametrize by the excess number of particles in the 
cluster, 



(4) 



where the cumulative mass is 



AN= / ipir)-p^)dr 



and, as pointed out previously, these need not be the 
same since one could increase AN while holding the ra- 
dius constant and increasing the width. This ambiguity 
is the fundamental problem with methods based solely on 
free energy considerations: there is no obvious method to 
determine which approach is preferable (or more perti- 
nently, which is chosen by Nature). 

Previously, it was shown that under the assumptions 
listed above, the MLP can be determined by gradient 
descent on the free energy surface. This means that first 
the saddle point is located and then one solves 



dxi 



(5) 



where the sign is chosen according to whether one is mov- 
ing uphill from the initial state, pi, to the saddle point 
(plus sign) or downhill from the saddle point to the final 
state, (minus sign). Alternatively, one can start at the 
saddle point and solve this equation with the minus sign 
to determine the two halves of the path. The informa- 
tion about dynamics is contained in the matrix of kinetic 
coefficients, g~j^, which are calculated as the inverse of 
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dm (r; x) dm (r; x) 



4ir2p(r;x) dxj 



dxj 



dr (6) 



m (r; x) = 47r / p{r )r dr . 
Jo 



(7) 



Another interpretation of Eq. ([5]) is that the MLP is de- 
termined by steepest descent on the free energy surface 
with the matrix gij playing the role of a Riemannian 
metric. In general, the probability density for any path, 
p (r; X (i))for < i < T, is given by 



P[x] =A/-exp( -^^[x] 



where the action is 



S[x]= [ Cix,±)dt, 
Jo 
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the Lagrangian is 
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dxi 
dt 
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(10) 

and where the normalization constant. A/", is independent 
of the path but otherwise unknown. Given two paths, 
X it) and y (t) , their relative probability can be calculated 
using these expressions provided they are parametrized 
in the same way. Note that Eq.® gives the probability 
density and not the probability of the path. The lat- 
ter would actually be P [x] ZJx where Z3x is the path 
measure. For example, if time were discretized using 
M + \ values ti = i (T/M) and if the parameters are 
then given by {tj) then the measure would be 

N M 

Dn. = H n dxij. This factor is irrelevant for comput- 

ing the relative probability of two paths described by the 
same parametrization since it would drop out of the ratio 
leaving the ratio of the probability densities, P [x] /P [y]. 
On the other hand, it is clear that one cannot meaning- 
fully compare two paths based on different parametriza- 
tion schemes since the measure, and the normalization 
factor, would in general be different. 

It may still be possible to compare the paths in an 
approximate manner by approximately translating one 
parametrization scheme into another. For example, a 
path based on a discretization of the density on a set of 
N lattice points, r^, cannot be directly compared to one 
based on a discretization over 2N lattice points, r[. How- 
ever, one can translate either parametrization into the 
language of the other by, e.g., using cubic-spline interpo- 
lation over the iV-point profiles to evaluate the density 
at the lattice positions of the 27V-point discretization. 

Another question is whether there is any way to com- 
pare the probabilities of pathways that are not deter- 
mined from the dynamical model? In particular, one 
would like to be able to evaluate the utility of the large 
number of existing heuristic approaches for the determi- 
nation of nucleation pathways which typically give a path 
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in terms of a reaction coordinate, which will be called s, 
rather than the time? For example, if the path is deter- 
mined by minimizing the free energy at constant excess 
particle number, N, for different values of N ranging from 
zero to the mass of the critical cluster, then the reaction 
coordinate, s, would be TV. It would seem that the only 
physical meaning of such a path is that a real system 
would follow it in the course of time so that if the path 
X (s) goes from the initial state at s = to the saddle 
point at s = Smax then, this must mean that the physical 



system begins at p {r;t = 0) = p (r; x (s = 0)) and ends 
at the saddle point at p{r]t = T) = p (r; x (s = Smax))- 
Since both time and the reaction coordinate vary mono- 
tonically during this processes, it must be that they are 
related so that we can write 



dp ds dp _ . . dp 
dt dt ds ds 



(11) 



with v{s) = ds/dt being the speed along the path. Writ- 
ing substituting this into the expression for the path 
probability wc get 



P [p] = TVcxp 
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ds dx, 
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We can now maximize the probability with respect to speed the system advances along the proposed path. Di- 



variations in v{s) to find that 

dx^ dxi dn _i ari 2 



rect evaluation of the second functional derivative of the 
path probability density shows that this is indeed a max- 
imum. The explicit form of the induced dynamics is 



or 
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This relation provides the desired expression for the The probability density itself becomes 
I 



P [p] = Mcxp 



dxi dxi 
ds ds 



9kj 



dx 



dxi dfl 
dxk ds dxi 



ds 
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which, with the time eliminated, gives a method to com- 
pare the probabilities for parametrized paths. Notice 
that this expression is invariant under a reparametriza- 
tion of the reaction coordinate, e.g. ds — >■ u{s')ds', so 
that it is a purely geometric quantity. In fact, if we fur- 
ther introduce the gradient force 
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dxi 



(17) 



it can be written as 



P [p] = A/'cxp 



dx 

ds 



Ibl sin 




where |b| = ^Jbigabj and where 0(^,b) is the angle 
between ^ and b. This form emphasizes the purely 
geometric nature of the optimized path. In particular, 
the parametrization-invariance means that the result is 



independent of the value of the upper limit, Smax, so that 
any convenient rescaling is allowed. 



III. APPLICATIONS 

The dynamical theory on which these results arc based 
is applicable to, e.g., globular proteins in solution which 
can be modeled at the crudest level as spherical molecules 
interacting via a short-ranged effective potential and sub- 
ject to Brownian forces due to the solvent. Here, the 
model potential of ten Wolde and Frenkel consisting of a 
hard core and short-ranged attraction. 
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FIG. 1. The spherically-symmetric concentration (density) 
distribution at various points along the classical path, as de- 
termined from Ea. (|20|) . The initial state has uniform concen- 
tration, p(r) — 0.075, and the value of the concentration at 
the origin increases monotonically along the nucleation path- 
way. 



will be used. The energy scale is set by e and the hard- 
core radius is a while the parameter a determines the 
distance of the attractive minimum from the hard core: 
the typical value for globular proteins of a = 50 will 
be used. The bulk free energy is approximated using 
thermodynamic perturbation theory. 

In CNT, it is assumed that the material in the in- 
terior of a cluster is in the bulk state: e.g., if a liquid 
phase is being nucleated from a gas then it is assumed 
that the density inside the cluster is always that of the 
bulk liquid and the cluster grows simply by increasing 
its radius from zero. More detailed non-classical results 
allow for the possibility of the interior density to vary 
from the bulk value and the typical behavior, as shown 
in Fig. 1, is for the density to begin at that of the back- 
ground gas and to increase as the radius increases until 
eventually reaching that of the bulk when the cluster is 
very large. This is associated with the fact that in con- 
trast to CNT the interface between the cluster and the 
gas is of finite extent so that when the cluster is small, 
all molecules may be considered to be in an intermedi- 
ate, interfacial region. Recently, it has been shown that 
solving the equations for the MLP yields an unexpected 
result: the cluster does not begin with a small, local- 
ized increase in density but rather it starts as a spatially 
extended density fluctuation with the actual nucleation 
event occuring within this structure. The initial radius is 
infinite and the excess mass is finite so that the density is 
that of the vapor. The first part of the process involves 
the gradual decrease in the radius with the excess mass 
remaining nearly constant so that the density increases 
slowly. This represents the formation of a density fluctu- 
ation containing excess mass relative to the background. 
The second stage of the process is the formation of a nu- 



FIG. 2. Some aspects of the nucleation pathway as deter- 
mined using a typical DFT method, left panel, and deter- 
mined from the dynamical theory, right panel. The figures 
show the central density, p(r = 0), excess particle number, A*', 
and excess free energy, Afi, as functions of distance along the 
nucleation pathway. The left panal, labeled "Classical", is 
the result of the heuristic theory, Ea. (|20|l . whereas the right 
panel shows the most likely path (MLP), as described in the 
text. All quantities, including the ordinate, have been scaled 
by their values at the critical cluster. (Note that despite the 
differences in the paths, the critical cluster is uniquely deter- 
mined by the free energy and so is the same for both paths.) 
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FIG. 3. The same as Fig. ([2]) but showing the equimolar ra- 
dius, Rcq and the "total" radius, i?totai, as a function of dis- 
tance along the nucleation pathway for the path determined 
from DFT gradient descent, left panel, and from the dynam- 
ical theory, right panel. 



clcus within this region of slightly enhanced density. The 
excess mass within the density fluctuation is the basis for 
the formation of the cluster which then goes on to grow 
as in the classical scenarios. 

In order to make a quantitative comparison of these 
different pathways, calculations of the most likely path 
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for the nucleation of high-concentration droplets in a 
low-concentration protein solution were performed as de- 
scribed in Refi and using a "classical" method which is 
known to give a good quantitative description of con- 
strained clustersSS. The determination of the MLP 
amounts to the solution of Eq.©. The classical calcu- 
lation takes a very similar form to the MLP calculation 
despite the fact that it was proposed on purely heuristic 
grounds prior to the MLP method: it involves gradient 
descent on the free energy surface 



dxi 
ds 



1 
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(20) 



but with a heuristic metric given by 



ffy (x) = 



dp (r; x) dp (r; x) 



dxi 



dxj 



dr. 



(21) 



Note that the first of these equations is equivalent to 
Eq.(l5]), but with distance along the path used as the 
independent variable rather than "time" (since, in the 
classical theories, there is no dynamics and so no natu- 
ral concept of time). The similarity between these two 
calculations is striking and one of the advantages of the 
dynamical approach is that it shows that gradient de- 
scent on the free energy surface does indeed characterize 
the nucleation pathway provided the correct metric is 
used. The "correct" metric is basically defined in terms 
of mass differences rather than density differences which 
can be traced to the fact that the underlying theory is 
based on fluctuating hydrodynamics in which mass is 
conserved. In contrast, the heuristic metric was simply 
guessed based on the criterion of simplicity and a preju- 
dice towards the use of local density as the fundamental 
variable and, so, no physical meaning can be attached to 
the paths derived from it. 

In all calculations, spherical symmetry was assumed so 
that the configuration at any given instant is character- 
ized by the variation of density as a function of distance 
from the origin. Figure [1] shows the density distribu- 
tion taken from a sequence of points along the classical 
pathway. The concentration begins as a constant, equal 
to the concentration of the solution. The formation of 
a cluster involves a monotonic increase in density near 
the origin until a density close to that of the bulk high- 
concentration solution is obtained. Beyond this point, 
the cluster grows via a monotonic increase in its radius. 

For the MLP it is known that the path probability 
(using either Eq.® or Eg. (1161) ) is simply given by 



MLP 



(22) 



where Af2 is the difference in free energy between the be- 
ginning and end points. Numerical evaluation of Eqs.®- 
(jl6p confirms this relation to a high degrees of numeri- 
cal accuracy. The "classical" model gives the sequence 
of profiles shown in Fig. 1. The log of the path prob- 
ability generated by this ansatz is shown in Figure 2. 
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FIG. 4. The action for the classical and MLP nucleation 
pathways shown as a function of excess free energy. Note 
that in both cases, the path probability density is related to 
the action by P ~ e~^. Also, both paths end at the critical 
cluster which is the same for all paths. The curve for the 
MLP is trivial as one knows^ that Smlp = —Afl. 



The energy barrier in these calculations is (3 Aft ^ 88 so 
that Pmlp — Me~^^ . In contrast, the probability of the 
heuristic path is Pdassicai — ■N'e~^^^ so that it is as if 
the energy barrier were 480fcsr rather than SSksT. The 
relative probability is 



Pcla 



ssical 



-398 n 

^MLP 



(23) 



so that the heuristic path is very unlikely compared to 
the MLP. This simple comparison illustrates the fact that 
the path probability is quite sensitive to the path and is 
not determined solely by the free energy barrier. 

Finally, the variation of the relative path probabili- 
ties as a function of the supersaturation is illustrated 
in Fig. [5] which shows the difference in the action 
as a function of size of the critical cluster. At least 
for the range of cluster-sizes shown in the figure, there 
is an almost linear relation which implies the relation 

Pclass^cal ^ C"^-"-' /^p^^^ ^ 



IV. CONCLUSIONS 

Based on the recently proposed dynamical approach to 
nucleation, a method has been derived for determining 
the maximum probability of observing any given nucle- 
ation pathway. The method is applicable to all pathways 
without regard to how they are constructed. This allows 
one to then ask whether hcuristically-derived pathways 
are reasonable approximations to the most likely path- 
way. The formalism has been illustrated for the problem 
of the nucleation of high-concentration droplets from a 
low-concentration solution of globular proteins. It was 
shown that the probability for a "classical" nucleation 
pathway, which begins with a small, localized cluster. 
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FIG. 5. The difference in action between the classical path 
and the MLP, AS = Sdassicai — Smlp, for different val- 
ues of supersaturation corresponding to different critical clus- 
ter sizes. The relative path probability for the two paths is 
Pciassicai/ Pmlp = e^^'^. A linear fit to the calculated points 
gives A5 « 0.34iVc. 

compared to that of the most likely path, which invokes 
a two-step mechanism whereby nucleation begins with a 
spatially-extended, small-amplitude density fluctuation, 
scales like e~^=''"""'/'^ so that for all but the smallest 
clusters, it is extremely small. This comparison serves to 
show that classical paths, even though intuitively appeal- 
ing, can have very low probability of occurrence. Con- 
trary to the expectation that all that matters for nu- 
cleation is the free energy barrier, these results show 
that there is a dramatic quantitative difference in the 
likelihood of observing the classical scenario compared 
to the non-classical one. Of course, only one family 
of "classical" pathways has been investigated here and 



there are many alternatives such as those found in Rcfs. 
[13, EH andm. It is possible that some of these com- 
pare more favorably with the non-classical Most Likely 
Path, but this can only be resolved by direct calculations. 



ACKNOWLEDGMENTS 

This work was partially supported in part by the Eu- 
ropean Space Agency under contract number ESA AO- 
2004-070 and by FNRS Belgium under contract C-Net 
NR/FVH 972. 

^J. D. Gunton, A. Shiryayev, and D. L. Pagan, 
Protein Condensation: Kinetic Pathways to Crystallization and Disease 
(Cambridge University Press, Cambridge, 2007). 

2p. G. Vekilov, Crys. Growth and Design 4, 671 (2004). 

^K. G. Libbrect, Rep. Prog. Phys. 68, 855 (2005). 

^J. Bernstein, Polymorphism in Molecular Crystals (Oxford Uni- 
versity Press, Oxford, 2002). 

'""T. Thrclfall, Organic Process Research and Development 7, 1017 
(2003). 

•^I. Stranski and D. Totomanow, Z. Phys. Chem. 163, 399 (1933). 
^L. O. Hedges and S. Whitelam, J. Chem. Phys. 000, 000 (2011). 
*J. F. Lutsko, J. Chem. Phys. 135, 161101 (2011). 
9j. F. Lutsko, J. Chem. Phys. 136, 034509 (2012). 
lOy. Talanquer and D. W. Oxtoby, J. Chem. Phys. 100, 5190 
(1994). 

"J. F. Lutsko, Europh ys. Lett. 83, 46007 (2008). 

12 J. F. Lutsko, J. Chem. Phys." 129, 244501 (2008). 

l^S. Ghosh and S. K. Ghosh, J. Chem. Phys. 134, 024502 (2011). 

^^M. J. Uline and D. S. Corti, Phys. Rev. Lett. 99, 076102 (2007). 

^^T. Philippe and D. Blavettc, J. Chem. Phys. 135, 134508 (2011). 

I'^A. J. Archer and R. Evans, Molecular Physics 109, 2711 (2012). 

i^D. Kashchicv, Nucleation : basic theory with applications 

(Butterworth-Heinemann, Oxford, 2000). 
i*C. Qiu and T. Qian, |J. Chem. Phys .' 131, 124708 (2009). 
i ^X. Cheng, L. Lin, W. E, P. Zhang, and A.-C. Shi, 

Phys. Rev. Lett? ] 104, 148301 (Apr 2010). 
^ffj. F. Lutsko, J. Chem. Phys. 134, 164501 (2011). 



